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I.  Introduction 


This  Annual  Scientific  Report  covers  work  performed  under 
Contract  No.  N00014-77-C-0102,  entitled  Excimer  Potential  Curves.  This 
report  describes  the  present  status  of  our  effort  to  develop  and  implement 
semi-empirical  and  theoretical  methods  for  obtaining  potential  curves  of 
diatomic  excimer  systems.  Our  emphasis  is  on  developing  and  testing 
methods  which  will  be  reasonably  accurate  yet  will  not  require  long  lead 
times  for  development  and  will  not  require  excessive  amounts  of  computer 
time  for  production  runs.  The  object  is  to  enable  experimentalists  to 
choose  or  reject  possible  laser  systems  on  the  basis  of  inexpensive 
theoretical  calculations  rather  than  on  the  basis  of  expensive  and  time- 
consuming  experiments. 

We  are  particularly  interested  in  developing  methods  that  are  applicable 
to  excimer  systems  because  of  the  current  emphasis  on  these  systems  as 
candidates  for  efficient,  high-power  visible  and  ultraviolet  lasers. 

After  consultation  with  A.  V.  Phelps  and  A.  Gallagher  of  J.I.L.A.,  we 
have  decided  to  concentrate  initially  on  molecules  of  rare  gases  with 
Thallium,  Indium  or  Gallium. 

For  systems  with  a large  number  of  electrons,  such  as  Tl-R  , present 
Configuration  Interaction  (Cl)  programs  are  inadequate  and  offer  no  hope 
of  extension  in  the  immediate  future.  Consequently,  our  effort  is  con- 
centrated on  developing  and  testing  semi-empirical  methods  that  can  easily 
and  rapidly  be  applied  to  the  excimer  systems  of  interest.  This  report 
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covers  our  progress  in  three  areas:  1)  We  have  performed  a configuration 
interaction  calculation  on  GaKr  , 2)  We  are  applying  the  Gordon-Kim 
theory  to  the  group  1 1 IB-rare  gas  systems,  3)  We  aie  continuing  the  de- 
velopment of  the  semi-empirical  electron  scattering  model. 

Our  work  in  the  first  area  is  close  to  completion.  The  Cl  calculation 
on  GaKr  was  done  so  that  we  can  compare  the  results  of  our  approximate 
methods  with  accurate  potential  curves.  GaKr  was  chosen  as  the  test 
system  since  it  is  the  largest  of  the  systems  we  are  interested  in  for 
which  a Cl  calculation  can  be  done.  The  GaKr  curves  have  been  extrapolated 
to  give  curves  for  InKr  and  TlKr  ; however,  these  extrapolation  are 
very  approximate.  In  particular,  the  positions  and  depths  of  the  potential 
wells  for  InKr  and  TlKr  are  not  given  accurately.  These  potential 
curves  have  been  used  to  predict  the  spontaneous  emission  and  absorption 
coefficients  for  GaKr  , InKr  and  TiKr  . The  details  of  this  work  are 
given  in  Section  II,  which  is  a rough  draft  of  a paper  on  this  work  that 
will  be  submitted  for  publication. 

Although  at  the  time  the  proposal  was  written,  applications  of  the 
Gordon-Kim  electron  gas  theory  to  open  shell  systems  had  failed  to  give 
reasonable  results,  recent  applications  of  the  G-K  theory  to  open  shell  - 
closed  shell  interactions  have  been  successful  (M.  J.  Clugston  and  R.  G. 
Gordon,  JCP  66^,  239  (1977)).  Consequently,  we  are  now  investigating  the 
possibility  of  using  this  method  to  obtain  aoproximate  potential  curves 
for  GaR  . Our  preliminary  calculations  on  GaKr  are  described  in 


Section  III. 
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The  basic  theory  and  computational  method  for  the  electron  scattering 
model  are  described  in  Section  IV.  We  are  proceeding  with  this  approach, 
in  addition  to  the  Gordon-Kim  theory,  because  the  electron-scattering  model 
should  provide  a more  accurate  description  of  atomic  distortion  and  should 
be  applicable  to  a wider  range  of  molecules  than  the  G-K  theory. 


^ ^ a 
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II.  Electronic  States  of  GaKr  : Ab  initio  calculations  of  a prototype 
for  TlKr 

The  work  in  this  section  was  performed  in  conjunction  with  Thom  H. 
Dunning,  Jr.  of  Los  Alamos  Scientific  Laboratory  and  is  a rough  draft  of 
a paper  that  will  be  submitted  for  publication. 


Electronic  States  of  GaKr  : Ab  initio  calculations  of  a prototype  for 
TlKr 

I.  Introduction 

Among  the  metal-rare  gas  eximers,  the  TlXe  system  is  believed  to 
be  an  excellent  candidate  for  an  efficient  visible,  high  power,  tunable 
laserJ  Since  ab  initio  calculations  on  this  system  are  beyond  the  scope 
of  present  cDmputer  programs,  we  present  here  a configuration  interaction 
(Cl)  calculation  on  GaKr  , which  is  the  largest  group  IIIB-rare  gas  system 
for  which  Cl  calculations  can  be  done.  We  use  the  calculated  GaKr  curves 
to  model  the  potential  curves  for  InKr  and  TlKr  . Although  our  model 
does  not  allow  further  extrapolation  from  TlKr  to  TlXe  , it  is  hoped 
that  these  calculations  will  yield  some  insight  into  the  properties  of  the 
TlXe  eximer.  In  addition,  Gallagher  has  recently  raised  the  possibility 
of  using  GaXe  as  a laser  if  the  6a  can  be  obtained  from  dissociation 
of  Gal.,  . Consequently  the  GaKr  curves  should  also  be  of  intrinsic 
interest. 

In  this  paper,  the  Cl  calculations  on  GaKr  are  presented  along  with 
the  model  calculations  on  InKr  and  TlKr  . These  potential  curves  are 
used  as  the  basis  of  a classical  calculation  of  the  emission  and  absorption 
coefficients  for  these  systems. 
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II.  Details  of  the  calculation 

A.  Basis  set 

The  calculations  use  (14sllp6d)  primitive  Gaussian  bases  for  gallium 
and  krypton"^  as  a starting  point.  The  core  orbitals  (Is ,2s ,2p,3s ,3p)  are 
sincly  contracted  to  the  Hartree-Fock.  atomic  orbitals  while  the  valence 
orbitals  (4s  and  4p)  and  the  3d  orbitals  are  each  described  by  two  contracted 
functions (see  Table  I).  The  resulting  (5s4p2d/5s4p2d)  contracted  bases 
are  constructed  using  the  general  contraction  scheme  of  Raffenetti.^ 

These  basis  sets  are  extended  to  include  polarisation  by  adding  two 
diffuse  s functions  (5s  and  5s')  and  a diffuse  p function  (5p)  to 
describe  the  lowest  Rydburg  orbitals.  The  exponents  for  these  orbitals 
(^(5$)  = .026,  c(5s')  = .011,  ^(5p)  =.01)  are  obtained  from  atomic  calcu- 
lations on  the  excited  states  of  Ga  . The  final  basis  set  thus  consists 
of  a (16sl2p6d)  primitive  basis  contracted  to  [7s5p2d]  for  Ga  and  a 
(14sllp6d)  primitive  basis  contracted  to  [5s4p2d]  for  Kr  . 

B.  SCF  calculation 

The  starting  point  for  the  Cl  calculation  is  a Hartree-Fock  calculation 
2 + 

on  the  T.  state 

....  13o^l40^77r^l5a^l60 


3 


The  inner  core  nolecular orbitals  (twelve  n , twelve  ti  and  four  fi) 
are  held  doubly  occupied  from  this  point  on  and  are  replaced  by  the 
rigorous  nonlocal  Hartree-Fock  potential. 

"core  = "n"  , 

i=core 

With  the  core  orbitals  removed  from  consideration,  it  is  convenient 
to  renumber  the  valence  orbitals  so  that  the  Hartree-Fock  configuration 
is  written  as 

1 2,  2.  2.  1 4 

1()  2a  3o  4alTT 

At  large  R the  correspondence  for  the  valence  orbitals  is 


la 

-> 

4sKr 

2a 

-> 

4pKr 

3a 

->- 

4sGa 

4a 

4pGa 

1 TT 

->■ 

4pKr 
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In  addition  to  the  valence  orbitals,  nine  o , six  tt  and  four  6 
virtual  orbitals  are  used  in  the  Cl  calculations.  The  lowest  virtual  orbitals 
(5o,  6a,  7a,  2v,  3tt),  which  correspond  to  the  Ga  5s,  5s',  4p  and  5p 
atomic  orbitals  for  large  internuclear  separations,  are  obtained  by  the 
improved  virtual  orbital  (IVO)  procedure.^  The  IVO  orbitals  are  obtained 
by  removing  the  electron  from  the  4a  (valence)  orbital  of  the  above 
configuration  and  calculating  the  virtual  orbitals  for  the  (N-1 )-electron 
Hami 1 tonian. 


C.  Cl  calculations 

Full  polarization  Cl  (POL-CI)  calculations®  which  provide  a balanced 
description  of  all  states  of  the  4s-4p  manifold  were  carried  out.  A set 
of  reference  configurations  was  chosen  (see  Table  II)  to  describe  the  dominant 
configurations  for  the  and  states  of  the  molecule  and  the 

state  of  the  ion. 


The  full  POL-CI  c<il culations  include  all  (1+2)  electron  excitations 
relative  to  each  reference  configuration  subject  to  the  restrictions  that 
no  more  than  one  electron  occupy  the  Rydberg  5o-  orbital  and  no  more  than 
one  electron  occupy  any  virtual  orbital  (6a,  2tt  etc.).  This  results  in 
764  spatial  and  2314  spin  configurations  for  the  ^e"^  states,  556  spatial 
and  1565  spin  configurations  for  the  states  and  368  spatial  and  558 
spin  configurations  for  the  ^e  state.  A total  of  15a  , 14t:  and  26 
occupied  and  virtual  orbitals  are  used  in  the  POL-CI  calculations. 
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III . Results  for  GaKr 

A.  Electronic  states  without  spin-orbit  coupling 

2 1 

The  interaction  of  ground  state  Ga  ( S)  and  Kr  ( S)  atoms  gives 
2 2 + 

rise  to  a n and  a z state  as  shown  schematically  in  Fig.  1,  In  the 
2 + 

I state  the  occupied  4p  Ga  orbital  is  oriented  along  the  molecular 

? 2 + 
axis  while  in  the  n state  it  is  perpendicular  to  the  axis.  The  Z 

state  should  be  essentially  repulsive  inside  the  van  der  Waals  region 

because  of  the  unfavorable  interaction  of  the  4pa  orbitals  on  Ga  and 
0 

Kr  The  n state  does  not  have  this  repulsive  interaction  so  it  should 

★ 2 +1 

be  more  attractive.  The  interaction  of  excited  Ga  ( S)  and  ionic  Ga  ( S) 

2 + 1 + 

atoms  with  ground  state  Kr  gives  rise  to  a z and  a z state 
respectively.  In  both  cases  these  states  are  expected  to  be  somewhat 


attractive. 
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The  potential  energy  curves  (Fi.j.  2 and  Table  III)  from  the  POL-CI 
calculations  generally  follow  the  behavior  predicted  from  these  theoretical 
considerations.  These  calculations  are  not,  however,  designed  to  treat 
long-range  dispersion  forces.  Many  excitations  which  contribute  to  a 
C^r  ^ attraction  are  not  included  in  the  wavef unction.  For  these  reasons, 
one  would  expect  deeper  wells  in  all  the  potential  curves  including  the 

two  lowest,  generally  repulsive,  states.  As  predicted, the  1 n state 

* 2 + 2 + 

IS  less  repulsive  than  the  1 Y.  state.  The  bound  excited  state  2 E 

has  a minima  at  6.36  a.u.  about  3.0  eV  above  the  ground  state  curve,  while 

the  minima  in  the  ion  curve  is  at  6.28  a.u.  and  5.6  eV  above  the  ground 

state.  The  well  dei  iis  are  compared  with  those  obtained  by  Gallagher^^ 

in  Table  IV,  and  as  expected,  Gallagher's  wells  are  deeper. 

The  dipole  moments  of  these  states  and  the  transitiuM  moments  between 

2 + 

z and  the  lower  states  have  also  been  calculated  and  are  given  in 


Table  V and  Figs.  3 and 
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3 Electronic  states  including  spin-orbit  coupling 

A complete  treatment  of  the  electronic  states  of  GaKr  must  include 
the  effects  of  spin-orbit  coupling.  The  states  considered  here,  which 
dissociate  to  the  closed-shell  ground  state  of  Kr  and  an  open-shell  state 
of  Ga  or  Ga^  are  influenced  only  by  •.  spin-orbit  matrix  e.ements  of 
the  open-shell  atom. 

Following  the  procedure  used  previously, we  have  adopted  a 

simple  model  for  including  the  effects  of  spin-orbit  coupling  on  the 

calculated  potential  energy  curves  and  wavefunctions . The  experimental 

spin-orbit  parameters  for  the  open-shell  atom  (Ga  and  later  In  , T1 ) are 

used  to  determine  the  matrix  elements  of  the  spin-orbit  interaction,  , 

coupling  the  molecular  states  at  infinite  separation  and  these  matrix 

elements  are  assumed  to  be  independent  of  internuclear  distance.  The 

resulting  spin-orbit  matrix  is  added  to  the  dianonal  matrix  of 

electronic  energies  H , (R)  = 6..E-(R)  : 

^ elec'  ' 1 j 1 

H(R)  = 

and  the  total  matrix  H is  then  diagonalized  at  each  internuclear  distance 
R . Thus,  in  addition  to  the  assumption  tnat  the  spin-orbit  matrix  elements 
do  not  change  as  functions  of  R , this  model  assumes  that  only  one-center 
terms  need  be  included  so  that  only  the  spin-orbit  coupling  on  Ga  is 
important.  We  would  expect  this  procedure  to  provide  reasonable  results 
to  the  extent  that  the  molecular  states  retain  the  identity  of  the  atomic 


4 


8 


states  from  which  they  are  formed.  The  spin-orbit  matrices  are  given 

in  Table  VI  along  with  the  atonic  parameters  used  in  these  calculations. 

2 2 

The  parameter  x is  chosen  so  that  the  atomic  ^2/2  ^1/2 

have  energies  of  +X  and  -2A  respectively. 

We  shall  label  the  molecular  states  using  the  convention  of  Hund's 
case  (c)  where  , the  projection  of  total  angular  momentum  along  the 
molecular  axis,  is  the  only  good  quantum  number,  n is  defined  as 

= A + , where  A and  are  the  orbital  and  spin  angular  momentum 

, 2 + 
projection,  respectively.  The  molecular  e states  have  only  a 

O 

fj  = 1/2  component,  while  the  n states  yield  a = 3/2  and  = 1/2 
state.  The  states  in  the  Q representation  are  labeled  according  to 
increasing  energy  by  a Roman  numeral.  So  the  Q = 1/2  states  are  designated 
as  11/2  ,111/2  and  the  n = 3/2  states  are  13/2,113/2  . 

The  coupled  states  are  expressed  as  follows: 

|I  l/2>  = Cj.|l^5:‘*’> 

ill  l/2>  = -CjlV>  + Cj.|l^n> 

1 1 3/2 > = |i^n> 

The  spin-orbit  coefficient  C^,  and  are  given  in  Table  VII.  The 
potential  curves  and  transition  moments  for  GaKr  including  spin-orbit 
coupling  are  given  in  Table  VIII  and  Figs.  5 and  6. 
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IV.  Extrapolations  to  InKr  and  TlKr 

The  potential  curves  for  InKr  and  TlKr  are  modeled  on  the 

GaKr  curves.  The  lowest  excitation  energies  and  the  ionization  potentials 

for  the  series  Ga  , In  , T1  are  given  in  Table  IX.  As  can  be  seen,  this 

series  does  not  form  a steady  progression.  In  has  a lower  ionization 

potential  and  lower  excitation  energies  than  Ga  , as  expected  for  a 

heavier  atom.  However,  T1  has  a higher  ionization  potential  and  higher 

excitation  energies.  This  is  due  in  part  to  the  presence  of  a filled 

12 

4f  shell  in  T1  and  the  larger  spin-orbit  effects.  These  effects 
should  be  considered  when  extrapolating  the  GaKr  curves  to  InKr  and 
TlKr. 


To  simulate  InKr  and  TlKr  , the  experimental  spin-orbit  parameters 
for  In  and  T1  are  used  to  couple  the  GaKr  curves.  The  curves  are 
also  shifted  to  give  the  correct  atomic  excitation  energies  at  R = « 

(see  Table  IX).  This  procedure  should  give  at  best  a qualitative  description 
of  states  of  InKr  and  TlKr  , since  the  non-spin-orbit  coupled  states 
are  expected  to  have  quantitatively  different  well  depths  and  equilibrium 
separations . 

The  effect  of  the  increasing  spin-orbit  perturbation  in  going  from 
Ga  to  In  and  T1  is  evident  in  the  calculated  curves  which  are  given 
in  Figs.  7 and  8.  Only  the  well  depths  and  positions  for  the  I 1/2  and 
II  1/2  states  are  affected  by  the  spin-orbit  coupling.  The  other  states 
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are  the  same  as  those  for  GaKr  except  that  they  have  the  correct 
asymptotic  spacing. 

The  mixing  parameters  from  the  spin-orbit  coupling  calculations  for 
TlKr  are  also  used  as  the  coefficients  of  the  GaKr  wave  function  to 
estimate  the  transition  moments  for  this  system.  The  TlKr  transition 
moments  are  given  in  Table  X and  Fig.  9. 
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V.  Absorption  and  stimulated  emission  coefficients  for  possible  laser  transitions 

The  interest  in  tie  group  I I IB-rare  gas  systems  arises  from  the  possibility 
of  their  use  as  visible  laser  systems.  In  order  to  judge  their  usefulness  as 
lasers  it  is  convenient  to  calculate  the  absorption  and  stimulated 

emission  g^(l’)  coefficients.  Obtaining  quantum-mechanical  results  for  these 
quantities  would  require  a complex  calculation  which  would  be  inconsistent  with 
the  extrapolations  used  to  obtain  the  InKr  and  TlKr  curves.  Consequently, 
we  have  used  Gallagher's  analysis  [17],  which  is  based  on  the  classical 
Frank-Condon  principle. 

In  order  to  obtain  g^  and  , the  Cl  curves  are  first  fit  by  Morse 
potentials.  The  parameters  for  these  Morse  potentials  are  given  in  Table  XI. 

These  parameters  can  then  be  used  in  Gallagher's  equations,  along  with  the 
atomic  transition  rate,  to  obtain  absorption  and  stimulated  emission  coefficients 
for  pressure  and  excitation  conditions  of  interest  to  experimentalists.  We 
have  calculated  these  coefficients  for  two  different  types  of  conditions. 

The  high  temperature  results  correspond  to  the  case  where  the  concentration 
of  the  metal  is  obtained  from  the  vapor  pressure  of  the  metal  itself,  while 
the  low  temperature  results  correspond  to  obtaining  the  required  concentration 
of  the  metal  from  vaporization  of  MI^  (M  = Ga,  In,  or  Tl).  This  latter 
condition  has  been  suggested  by  Gallagher  as  a possible  means  of  obtaining 
high  concentrations  of  the  metal  at  low  temperatures.  In  both  cases  the 
densities  used  are  10^^/cm^  for  Kr  , 10^^/cm^  = 3[M  " 1.5[M  ^P3/2- 

and  2 x 10^^/cm^  = [M  • The  resulting  absorption  and  stimulated 
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m 
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emission  coefficients  for  GaKr  , InKr  and  TlKr  are  cjiven  in  Figures  10 
to  15. 


m 
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Table  I.  Gaussian  exponents  and  contraction  coefficients 


Exponents 


Contraction  coefficients 


Gallium  atom 


Is 

2s 

3s 

4s 

4s' 

457600. 

.000222 

-.000069 

. 000026 

-.000006 

0.0 

68470. 

.001732 

-.000535 

. 000205 

-.000048 

0.0 

15590. 

.008952 

-.002814 

.001070 

-.000247 

0.0 

4450. 

.035874 

-.011275 

.004337 

-.001007 

0.0 

1472. 

.114000 

-.038495 

.014707 

-.003399 

0-0 

541.3 

.274138 

-.100714 

.039748 

-.009279 

0.0 

214.8 

.414793 

-.211832 

.084475 

-.019587 

0.0 

88.81 

.275395 

-.175448 

.079654 

-.019104 

0 . 0 

27.18 

.029561 

.479840 

-.291821 

.072753 

0.0 

11.54 

-.006815 

.634145 

-.527118 

.134137 

0.0 

3.303 

.002253 

.069592 

.583707 

-.181778 

0 . 0 

1.334 

-.001017 

-.012299 

.674103 

-.358241 

0.0 

.1947 

.000251 

.002774 

.028077 

.615164 

0.0 

.07158 

0.0 

0.0 

0.0 

0.0 

1 .0 

5s 

5s' 

026 

1.0 

0.0 

on 

0.0 

1.0 

2p 

3p 

4p 

4p' 

5p 

3274. 

.001513 

-.000576 

. 000094 

0.0 

0.0 

765.4 

.013070 

- . 004981 

.000800 

0.0 

0.0 

241.6 

.067263 

-.026421 

.004337 

0.0 

0.0 

89.39 

.219542 

- . 089529 

.014443 

0.0 

0.0 

36.36 

.421107 

-.186734 

.031377 

0.0 

0.0 

15.60 

.376515 

-.144494 

.021501 

0.0 

0.0 

6.472 

.089425 

.258956 

-.046233 

0.0 

0.0 

2.748 

-.000502 

.570187 

-.125293 

0.0 

0.0 

1.090 

.001761 

.325305 

-.045636 

0.0 

0.0 

.2202 

-.000247 

.016563 

.452811 

0.0 

0. 0 

.06130 

0.0 

0.0 

n.o 

1.0 

o!o 

.01 

0.0 

0.0 

0.0 

0.0 

1.0 

3d 

4d 

59.66 

.031949 

0.0 

17.10 

.163546 

0.0 

6.030 

.367457 

0.0 

2.171 

.456851 

0.0 

.6844 

.305161 

0.0 

.160 

0.0 

1.0 
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Table  I.  Gaussian  exponents  and  contraction  coefficients 


Exponents 

Contraction 

coefficients 

Krypton  atom 

Is 

2s 

3s 

4s 

4s' 

605700. 

.000231 

-.000073 

.000029 

- . 000009 

0.0 

90300. 

.001755 

-.000551 

.000221 

-.000070 

0.0 

20920. 

.009076 

-.002894 

.001159 

-.000369 

0.0 

5889. 

.036990 

-.011834 

.004781 

-.001522 

0.0 

1950. 

.116154 

-.039826 

.016056 

-.005118 

0.0 

718.2 

.278401 

-.104801 

.043454 

- . 01 3886 

0.0 

285.4 

.415746 

-.217093 

.091899 

-.029537 

0.0 

118.6 

.267204 

-.175562 

.083789 

-.027309 

0.0 

38.16 

.027870 

.471395 

-.303023 

.103498 

0.0 

16.45 

-.005998 

. 636794 

-.570620 

.208810 

0.0 

5.211 

.002217 

.082255 

.501751 

-.235737 

0.0 

2.291 

-.001092 

-.014138 

.760483 

-.553570 

c.o 

.4837 

.000306 

. 003289 

.044857 

.701123 

0.0 

.1855 

0.0 

0.0 

0.0 

0.0 

1.0 

2p 

3p 

4p 

4p’ 

4678. 

.001392 

-.000569 

.000156 

0.0 

1120. 

.011666 

-.004777 

.001286 

0.0 

357.1 

. 060858 

-.025631 

.007059 

0.0 

131.4 

.2’n040 

-.092159 

.024990 

0.0 

52.86 

.42i000 

-.200936 

.056870 

0.0 

22.70 

.383515 

-.160784 

.040225 

0.0 

9.547 

.097383 

.267789 

-.084756 

0.0 

4.167 

-.001087 

.585908 

-.240291 

0.0 

1.811 

.002209 

.291397 

-.038636 

0.0 

.5337 

-.000509 

.015484 

.599154 

0.0 

.1654 

0.0 

0.0 

0.0 

1.0 

3d 

4d 

125.6 

.019168 

0.0 

33.31 

.125638 

0.0 

12.15 

. 366069 

0.0 

4.350 

.502482 

0.0 

1.494 

.264377 

0.0 

.35 

0.0 

1.0 
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Table  II.  Reference  Configurations 


states 


lCT^2CT^3o^4al  TT 
1a^2a^3a^5al7r 
1a^2a^3a^6al'iT 


4 

4 

4 


2 


n states 


1a^2c;^3a^l/2TT 

la^2a^3a^l/3TT 

State  (GaKr^) 

1 2,  2-  2,  4 
la  2a  3a  lir 


.J 
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Table  III. 

POL-CI  calculations  on  the  low-lying  states  of  GaKr 
and  the  ground  state  of  GaKr"^  . All  energies  are 
relative  to  -4674.  hartrees. 

R 

2^1:* 

I'l* 

oo 

-1.200042 

-1.095050 

-1.198917 

-0.992860 

15.00 

-1.200121 

-1.095025 

-1.198978 

-0.993044 

10.00 

-1.200290 

-1.094061 

-1 .199328 

-0.994018 

8.00 

-1.198650 

-1.094202 

-1.200204 

-0.996703 

7.00 

-1.193963 

-1.095454 

-1.200450 

-0.999505 

6.00 

-1.179381 

-1.096768 

-1.197718 

-1.002018 

5.00 

-1.135164 

-1.088552 

-1.179245 

-0.992868 

4.50 

-1.086883 

-1.064404 

-1.148099 

-0.967583 

4.00 

-1.018697 

-0.992381 

-1.079327 

-0.906827 

3.75 

-0.968656 

-0.933786 

-1.025394 

-0.856247 

' ’*4*  4 
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Table  IV.  Potential  well  depths 


Mole  - State 
cule 


Cl  Morse  Fit  Gallagher 

R AF(eV)  R AF(eV)  R AE(eV) 


GaKr  Iy2i^y2) 

-v-.OlO 

7.55 

.021 

^3/2^h/2'> 

'\-.040 

7.18 

.040 

IIi/2(Ai/2) 

'x.lO 

.006 

10.15 

.00642 

6.36 

.064 

6.26 

.080 

GaKr 

6.28 

.26 

5.95 

.252 

I nKr  ^1^2 

'V.013 

7.86 

.00778 

I3/2' 

^^1/2 

MO 

^^.006 

9.98 

.0064 

IIIl/2^ 

TlKr  1^/2 

'\;8 

-v^.OlO 

8.5 

.012 

7.01 

.024 

^3/2^ 

6.58 

.062 

^^1/2 

MO 

'\<.006 

9.88 

.064 

IJIl/2^ 

lo^ 

6.09 

.107 

^same  as  GaKr 

Table  V. 

Dipole  and  transition  moments  for  the  low-lying  states  oT 

GaKr 

R 

x^n  2 V 2^5:’^-x^n 

15.00 

0.00598 

0.00156 

0.07558 

0.00183 

-1.29167 

-1.31011 

10.00 

0.02605 

0.02048 

0.56872 

0.01631 

-1.28976 

-1  29762 

8.00 

0.09152 

0.14168 

0.89914 

0.04761 

-1 .27767 

-1 .23098 

7.00 

0.18861 

0.30443 

0.95672 

0.07894 

-1.26770 

-1.19024 

6.00 

0.41543 

0.62092 

0.81671 

0.11663 

-1.25536 

-1.18926 

5.00 

0.90382 

1.20144 

0.45889 

0.09916 

-1.23487 

-1.44020 

4.50 

1.24437 

n qcojii 

W • ^ wv>  -r  t 

1.19373 

-0.13653 

-1.20897 

-2.06755 

4.00 

1.64453 

-1 .27575 

5.47300 

-0.84619 

-0.86801 

- .30554 

3.75 

1.88349 

-1.13125 

6.45764 

-0.90367 

-0.78769 

.40946 
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Table  VI. 

Quantities 

for  spin-orbit  matrices 

Q = 1/2 

Q = 3/2 

0 

/Ta 

2„ 

n 

x/Fa 

-A 

Ga  A = .001255  au 
In  A = .00336 
T1  A = .311835 
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Table  VII. 

Spin-orbit 

coefficients  for  the 

n = 1/2  states 

R(a„) 

4 

GaKr 

3.75000 

.99955 

.02998 

4.00000 

.99960 

.02813 

4.50000 

.99961 

.02787 

5.00000 

.99927 

.03811 

6.00000 

.99641 

.08468 

7.00000 

.98199 

.18892 

8.00000 

.93379 

.35782 

10.00000 

.82915 

.55902 

15.00000 

.81649 

.57736 

InKr 

3.75000 

.99704 

.07690 

4.00000 

.99738 

.07238 

4. 50 ',00 

.99742 

.07175 

5.0UU00 

.99534 

.09644 

6.00000 

.98063 

.19588 

7.00000 

.93712 

.34900 

8.00000 

.87677 

.48091 

10.00000 

.82132 

.57047 

15.00000 

.81649 

.57736 

TlKr 

3.75000 

.97505 

.22197 

4.00000 

.97731 

.21180 

4.50000 

.97763 

.21034 

5.00000 

.96502 

.26219 

6.000C0 

.91737 

.39803 

7.00000 

.86679 

.49867 

8.00000 

.83612 

.54854 

10.00000 

.81788 

.57540 

15.00000 

.81649 

.57736 

Table  VIII.  Rydberg- valence  transition  moments  in  GaKr 
(with  spin-orbit  corrections) 


R 

III  1/2 

Z 

- I 1/2 

(X,Y) 

III  1/2 

Z 

- II  1/2 

(X,Y) 

III  1/2  - I 3/2 

(X.Y) 

15.00 

-0.75641 

0.74574 

-1.06969 

-0.52733 

-0.91335 

10.00 

-0.72540 

0.75618 

-1.0759? 

-0.50983 

-0.91200 

8.00 

-0.44047 

0.84363 

-1.14948 

-0.32327 

-0.90345 

7.00 

-0.22486 

0.88026 

-1.16880 

-0.16935 

-0.89640 

6.00 

-0.10071 

0.88449 

-1.18499 

-0.07517 

-0.88767 

5.00 

-0.05489 

0.87255 

-1.43915 

-0.03328 

-0.87318 

4.50 

-0.05762 

0.85454 

-2.06674 

-0.02383 

-0.85487 

4.00 

-0.00859 

0.61353 

-0.30542 

-0.01727 

-0.61378 
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Table  IX.  Atomic  states  of  Ga  , In  , and  T1 


Excitation  Energy 

State 

Ca 

In 

T1 

cm"^ 

eV 

cm”^ 

eV 

cm 

eV 

1/2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

^3/2 

826.24 

.10241 

2212.56 

.274228 

7792.7 

.965840 

^1/2 

24788.58 

3.07234 

24372.87 

3.020814 

26477.5 

3.281665 

I.P. 

48380. 

5.9963 

46669.93 

5.784348 

49264.2 

6.105886 

Table  X.  Rydberg-to-valence  transition  moments  in  model  TlKr 
(with  spin-orbit  corrections) 


III  1/2  - 11/2  III  1/2  -II  1/2  III  1/2  - 13/2 


R 

Z 

(X,Y) 

Z 

(X.Y) 

(X.Y) 

15.00 

-0.75641 

0.74574 

-1.06969 

-0.52733 

-0.91335 

10.00 

-0.74665 

0.74590 

-1.06130 

-0.52476 

-0.91200 

8.00 

-0.67524 

0.75539 

-1.02925 

-0.49558 

-0.90345 

7.00 

-0.59354 

0.77699 

-1.03169 

-0.44701 

-0.89640 

6.00 

-0.47336 

0.81433 

-1.09099 

-0.35332 

-0.88767 

5.00 

-0.37761 

0.84264 

-1.38982 

-0.22894 

-0.87318 

4.50 

-0.43489 

0.83575 

-2.02130 

-0.17981 

-0,85487 

4.00 

-0.06471 

0.59985 

-0.29861 

-0.13000 

-0.61378 
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Table  XI. 

Morse  fitting  parameters: 
V(R)  ' V(»>)  = Dg[u^  - 2u] 

where  u 

= exp  [B(Rg  - R)] 

Molecule 

State 

“e 

R (a„) 

e 0' 

B(V') 

GaKr 

^1/2 

.000772 

7.56 

.724 

^3/2 

.00147 

7.19 

.722 

^h/2 

.000625 

10.14 

.471 

III1/2 

.00295 

6.26 

.853 

InKr 

h/2 

hn 

. 000286 

7.86 

.817 

”1/2 

IIIl/2^ 

. 000474 

10.24 

.483 

TlKr 

h/2 

h/2 

. 00439 

8.52 

.633 

A 

^h/2 

III^/2^ 

.000408 

10.12 

.502 

^Same  as  for  GaKr 


Fig.  1.  Orbital  diagrams  for  the  electronic  states  of  Ga  + Kr  and  Ga  + Kr. 

The  lobes  and  circles  represent  the  in-plane  and  out-of-plane  p orbitals 
the  dashed  circle  denotes  the  Rydberg  orbital. 
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Figure  4. 


DIPOLE  TRANSITION  MOMENTS  AMONG 
THE  LOW-LYING  STATES  OF  GaKr 
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Figure  5. 

THE  LOW-LYING  STATES  OF  GoKr  AND 
GoKr'*'  WITH  SPIN-ORBIT  CORRECTIONS 


I 


I 


1.5  2.0  2.5  3.0  3.5  4.0  4.5  5.0  5.5  6.0 
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F ' Figure  6. 
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DIPOLE  TRANSITION  MOMENTS  AMONG 
THE  LOW-LYING  STATES  OF  GoKr  WITH 

SPIN-ORBIT  CORRECTIONS 

I 


I 


2 


0 

-I 


Q 


32 


Figure  7. 


I 

k 


THE  LOW-LYING  STATES  OF  InKr  AND 
InKr"^  WITH  SPIN- ORBIT  CORRECTIONS 
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I Figu>^e  8. 


‘ THE  LOW-LYING  STATES  OF  TIKr  AND 

f TIKr^  WITH  SPIN* ORBIT  CORRECTIONS 
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Figiire  9. 


DIPOLE  TRANSITION  MOMENTS  CONNECTING 
THE  LOW-LYING  STATES  OF  TIKr  WITH 
SPIN-ORBIT  CORRECTIONS 


550  650 
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III.  Application  of  the  Gordon-Kim  theoiy  to  the  group  IIIB-rare  gas  systems 

In  recent  years,  the  electron  gas  methods  developed  by  Gaydaenko 
and  Nikulin  [1],  and  Gordon  and  Kim  [2]  and  modified  by  Rae  [3]  and  Cohen 
and  Pack  [4]  have  proved  successful  in  calculating  the  interaction  energies 
of  pairs  of  closed  shell  atoms  or  ions  [1-5],  of  atom-molecule  systems[6-7] , 
of  molecule-molecule  systems  [8],  and,  recently,  of  systems  with  one  open 
shell  atom  [9].  These  interaction  energies  are  much  more  accurate  than  the 
simplicity  of  the  Gordon-Kim  (GK)  method  would  lead  one  to  expect. 

The  recent  success  of  the  GK  method  in  calculating  closed  shell  - open 
shell  interaction  energies  has  prompted  us  to  use  this  method  to  calculate 
the  interaction  of  group  1 1 IB-noble  gas  pairs.  The  GK  theory  has  been 
most  successful  in  cases  where  the  interaction  is  non-covalent  and  where 
the  atoms  are  relatively  undistorted.  For  these  reasons,  systems  such  as 
GaKr  would  seem  to  be  ideal  candidates  for  a GK  calculation. 

The  electron-gas  theory  is  briefly  reviewed  in  part  A,  and  our 
preliminary  results  on  GaKr  are  given  in  part  B. 

A.  The  electron  gas  method 

The  method  we  have  used  is  that  of  Gordon  and  Kim  [2]  as  modified  by 
Cohen  and  Pack  [4].  A more  detailed  description  of  the  theory  is  available 
in  these  two  papers.  Briefly  the  GK  theory  approximates  the  intermol ecnlar 
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potential  V(R)  at  the  distance  R by 
V(R)  V^^R)  = V^^(R)  + 


where 


''Sf  = 

and  these  three  terms  represent  the  kinetic,  Coulomb  and  exchange  inter- 
action energies,  respectively.  To  calculate  these  interaction  energies, 
the  electronic  charge  density  p is  approximated  by  the  sum  of  the  atomic 
charge  densities, 

p =.  p,  * ''b  • 

With  this  approximation,  the  Coulomb  interaction  can  be  calculated  directly, 
but  the  other  terms  are  all  estimated  by  the  formulas  for  the  energy  density 
of  a uniform  electron  gas  [4].  An  additional  modification  [3,4]  is  made 
to  the  exchange  energy  to  avoid  self-exchange  contributions. 

B.  Calculations 

We  have  modified  the  molecule-molecule  GK  interaction  program  of 
Parker,  Snow  and  Pack  [8],  which  allows  for  non-spherical  potentials,  to 
calculate  interaction  energies  for  closed  shell  - open  shell  atomic  pairs. 
The  density  of  the  open  shell  atom,  in  this  case  Ga  , is  divided  into 
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the  spherically  symmetric  core  density  plus  the  valence  density.  For  Ga 
in  the  ground  state  (4s  4p  ),  the  valence  density  Is  constrained  to  be  in 
a p orbital  directed  along  the  internuclear  axis,  to  form  a E molecular 
state,  or  perpendicular  to  it  to  form  a n state. 

In  order  to  expedite  the  evaluation  of  integrals,  the  atomic  density 
is  expanded  in  a set  of  Slater  basis  functions  with  the  coefficients  determined 
by  a least-squares  fit  [4].  The  basis  set  expansions  for  the  density  are 
then  read  into  the  GK  program  which  calculates  the  interaction  energy 
by  three-dimensional  numerical  quadriture. 

We  have  obtained  numerically  tabulated,  relativistic  Hartree-Fock 

densities  for  the  group  IIIB  and  rare  gas  atoms  from  Joseph  Mann  [10]. 

2 1 2 1 

The  densities  of  the  4s  4p  and  4s  5s  states  of  Ga  and  the  ground  state 
of  Kr  were  fit  with  small  sets  of  Slater  functions.  At  this  time  the 
fit  of  the  basis  set  expansion  is  not  very  good  ('v  10%). 


The  basis  set  expansions  obtained  from  this  fitting  procedure  were 
then  used  to  calculate  the  GK  interaction  energies.  Because  the  fitting 
procedure  does  not  normalize  the  density,  the  inaccuracy  of  the  present 
expansions  results  in  a spurious  Coulomb  repulsion  between  the  two  atoms 
We  hope  that  improving  the  basis  set  fit  will  correct  this  error. 
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The  GK  potential  curves  are  compared  with  the  Cl  results  (see  Section  II) 
in  Figures  1-3.  The  energies  are  plotted  with  respect  to  the  asymptotic 
energy  of  each  state.  The  curve  labeled  GKR  includes  Rae's  correction  while 
the  GK  curve  is  the  unmodified  GK  theory.  The  GKR  curve  for  the  state 
agrees  remarkably  well  with  the  Cl  result.  Unfortunately,  the  agreement  for 
both  the  1 n and  2 E states  are  not  as  good.  Specifically  the  depths  and 
positions  of  the  potential  wells  are  not  predicted  accurately.  We  are 
currently  working  on  improving  the  basi:,  set  fit  and  thereby  the  interaction 
energies. 
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IV.  The  Electron  Scattering  Model 

A.  Introduction 

This  section  describes  the  theoretical  basis  of  our  calculations 
using  the  electron  scattering  model.  For  an  excimer  system  AB  , where 
A is  a closed-shell  system,  most  of  the  states  of  interest  correspond 
to  the  asymptotic  situation  where  B is  excited  but  A is  in  its  ground 
state.  Fundamentally,  what  the  electron  scattering  theory  says  is  that  any 
charged  particle  in  B sees  a potential  , due  to  the  closed-shell 

system  A , which  is  the  same  as  if  the  charged  particle  were  scattered 
off  of  A . This  scattering  potential  is  corrected  for  the  fact  that  A 
is  "de-polarised"  relative  to  the  scattering  problem,  due  to  the  presence 
of  the  nucleus  and  other  electrons  of  B . This  theory  is  based  on  the 
model  interaction  potentials  and  response  functions  that  arise  out  of  the 
many  body  theory  (using  Schwinger  Functional  Derivatives)  and  that  have 
been  applied  to  scattering  problems  [1]. 

In  part  B the  basic  equations  are  presented  and  pos.  ble  methods 
for  calculating  the  interaction  energy  are  discussed  in  part  C.  The  use 
of  semiempirical  forms  for  the  effective  notentials  are  discussed  in 
part  D. 


The  Basic  Equations 
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Using  many-body  field  theoretic  methods  it  has  been  shown  that  the 
change  in  energy,  e , resulting  from  the  addition  of  an  electron  to  a 
closed-shell  reference  system  (referred  to  here  as  A)  is  given  by  the  one- 
particle  Dyson  equation 


T(r)i|)(r)  + Jdr'  ?:^(r;r ' :c)4>(r)  = c4(r) 


(1) 


where,.  T is  the  kinetic  energy  operator  and  (f)  is  the  Dyson  amplitude 
with  r and  r'  being  space-spin  coordinates.  Thus,  the  problem  reduces 


to  an  effective  one-particle  problem  in  which  this  particle  experiences 


an  effective  potential ^ E , which  represents  all  the  other  particles 
collectively,  taking  into  account  all  effects  such  as  polarization. 


correlation  and  exchange,  etc.  As  might  be  expected,  the  cost  of  this 

,A  . 


simple  formulation  is  that  r.  is  an  extremely  complicated  entity  which 
is  both  nonlocal  and  energy  dependent  and  which  cannot  rigorously  be 
brought  into  closed  form.  However,  it  has  been  possible  to  develop 


iTiiii  if  iiiaiiliVnii'lii  - 
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excellent  closed  form  approximations  to  this  potential  which  are  based  on 
well  founded  physical  concepts.  Most  notable  among  these  is  the  Random 
Phase  Approximation  (RPA)  potential,  , which  has  been  very  successfully 

usv?d  ; calculating  the  ionization  potentials,  excitation  energies, 
oscillator  strengths,  and  elastic  -scattering  phase  shifts  for  He  [21. 

Zpp^  has  also  been  used  to  accurately  calculate  the  ionization  and  ex- 
citation energies  of  Li  [3].  Moreover,  it  has  been  shown  that  this  ab 
initio  potential  encompasses  other  phenomenologically  derived  semiempirical 
potentials  which  have  bt  n used  by  other  workers  with  great  success  [4]. 

It  has  also  been  recently  demonstrated  that  by  applying  the  same  many- 
body  techniques  to  the  problem  of  two  electrons  added  to  a closed-shell 
reference  system  (A)  one  obtains  an  effective  two-particle  equation  for 
the  resulting  change  in  energy  [5] 

2 

.'"2)  , (2) 

where  in  addition  to  the  individual  one-particle  potentials  there  now 

A 

appears  a two-partic''c  effective  potential,  W , which  represents  how  the 

presence  of  one  particK  affects  the  potential  seen  1 the  other  particle 

and  which  reflects  the  fact  that  system  A can  act  as  a dielectric  medium 

to  shield  the  Coulombic  interaction  between  two  charged  particles.  As 
A 

with  z , this  potential  cannot  be  rigorously  brought  into  closed  form. 
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Hovjever,  an  excellent  ^ initio  closed  form  approximation  to  this  potential 
can  be  obtained  with  the  Random  Phase  Approximation,  , which  is  completely 

compatible  with  the  similarly  obtained  one-particle  potential,  . 


If  we  now  proceed  to  the  case  of  adding  m electrons  to  our  closed- 
shell  reference  system  it  follows  by  induction  and  can  be  proven  formally 
that  the  change  in  energy  is  given  by 


dr:  dr:  dr'  u''(r,  ,r . ,r^;r:  ,r:  ,r' ir)  + ...  +/dr;  ...  dr;v''(r, ..  ,r^;r] ..  .r;:^] 


where  our  notation  for  three-particle  and  higher  potentials  is  obvious. 

In  the  above  we  only  considered  the  addition  of  electrons  to  our  closed- 
shell  reference  system  A.  However,  we  can  also  add  nuclei  as  well;  the  only 
difference  being  that  when  acting  on  a nuclear  coordinate  all  potential 
terms  involving  electron  exchange  must  of  course  vanish.  Since  this  can 
be  trivially  accomplished  by  choosing  an  appropriate  (’  ificial)  nuclear 
spin  coordinate,  we  can  immediately  generalize  (3)  to  include  both  nuclei 
and  electrons  by  allowing  the  particles  to  carry  different  charges,  z , 

(where  z = -1  for  an  electron)  to  give 
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[.l;{T(r,)  - z.  /dr;  /(r.irjre))  t ^|z,z.(|-l-|  + /dr:dr/(r,  ,rj  irl  .rj 


)) 


(4) 


We  of  course  cannot  solve  (4)  since  the  potentials  involved  cannot  be  written 
in  closed  form.  However,  we  can  replace  these  potentials  by  their  RPA 
approximates  which  are  in  closed  form.  Furthermore,  realizing  that  we  are 
deriving  a theory  for  intermolecular  forces  which  is  essentially  perturbative 
in  nature,  we  will  now  assume  that  all  three-particle  and  higher  potentials 
can  be  neglected  so  that  our  equa,tion  becomes  simply 


= • (5) 


Since  our  interest  lies  in  the  calculation  of  potential  energy  surfaces, 
what  we  really  want  is  an  equation  for  the  change  in  energy  when  m^ 
electrons  and  m|^  nuclei  are  added  to  A with  these  nuclei  held  fixed 
at  specific  points  in  space  (which  we  will  refer  to  collectively  as  R 
with  the  actual  spatial  coordinate  of  nuclei  i being  R^  ).  That  is, 
we  want  to  be  able  to  separate  electronic  motion  from  nuclear  motion  so 
that  this  quantity, <£'(R)  , will  be  the  total  potential  experienced  by 
these  nuclei  and  will  satisfy  the  equation 

'"n 

lEt(R.)  +<f(R)]  X(R)  = eX(R)  (6) 

i--l  ^ 
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where  X(R)  is  a function  of  the  nuclear  coordinates  only.  As  it 
stands,  (5)  precludes  such  a separation  because  o'  the  energy  dependence 
of  the  potentials  involved.  To  overcome  this  problem  we  will  now  assume 
that  for  those  solutions  to  (5)  we  seek  this  energy  dependence  is  not 
strong  Furthermore,  we  will  also  assume  that  the  response  of  A to  the 
added  particles  is  instantaneous.  In  this  way  we  can  replace  the  non- 
adiabatic  energy  dependent  potentials  in  (5)  with  their  hermitian  energy 
independent  adiabatic  approximates  such  as  those  given  in  ref.  5.  With 
these  substitutions  (5)  does  become  separable,  and  by  taking 


- EXj(R)T.(r,...r^  ) ; <».|Tj 


> = 6.  . 

U 


(7) 


I S' 


Where  'i’(ri...r^  ) is  a function  of  only  the  electronic  coordinates 
we  find  that  <1?(R)  is  given  by 


= e^(r)  -_Ez.  j:^(r.)  + j:_z.z.(tr^4^+  W^(R.,R.))  , 


a 

where  E is  given  by  the  equation 


(8) 


P "n 

[H(ri...r  ;R)  + E ( fdr'.Z^Crjr' ) - Vz.  fdr:w^(r.  ,R.  jr! )) 

I nig  j 1 J J 1 1 J 1 

+ E ]'dr:dr^W^(r.,rj;rlrj}]  'l'(r^...r^  ) = E^(R)  'l'(r.,...r^  ) (9a) 


where  H is  the  usual  hamiltonian  for  m^  electrons  in  the  field 
of  m^  fixed  nuclei 
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;R) 


. (T(r.)  - 1 ^ E 

1=1  j=l  ' 1 j'  i>j 

"’e  -"e  JL 

Eh(r.)  + E 

1=1  ^ i>j  ^ 


r^-rjl 


(9b) 


and  where  all  one-  and  two-particle  effective  potentials  are  taken  to  be 
hermitian  adiabatic  approximates  to  the  true  potentials.  Note  also  that 
all  potential  terms  involving  the  m^  nuclei  are  now  written  as  local 
quantities  thereby  taking  into  account  the  previously  mentioned  fact 
that  there  are  no  exchange  (i.e.  nonlocal)  terms  in  the  potentials  when 
nuclei  are  being  considered. 


If  we  now  collectively  refer  to  these  m^  electrons  and  m^ 
nuclei  which  have  a fixed  internuclear  geometry  as  being  system  B , 
then  £‘{R)  is  the  change  in  energy  resulting  from  the  creation  of  system 
B in  the  vicinity  (as  measured  by  R ) of  A . If  R ” then 

which  is  just  the  energy  B itself.  Therefore,  the  intermolecular 
^ential  of  the  system  A-B  as  a function  of  the  separation  between 
A and  B is 

V*(R)  . E*(R)  - E°  . (10a) 


where 


V 


A 

B.nuc 


(R) 


(10b) 
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and  where 


= + 

^ B 


i>J 


i J |R^ 


■'ji 


(lOc) 


is  just  the  electronic  energy  of  isolated  system  B and  where  Eg(R) 
is  given  by 


[H(r^...r  ;R)  + V /dr  I (z^(r . ;r '. ) - f zy  (r . ,R.  ;rl ) ) 

i=l  ^ ^ ^ j = l J T J 1 


^e 

+ E /dr:drjW^(r^.  ,rj;rjrj)]  '•'gCr^.-.r^  ) = Eg(R)  '1'^  ]•  • • ) (lOd) 


Therefore,  given  that  we  know  Eg  and  that  we  have  available  good 
closed  form  adiabatic  approximates  for  and  our  problem  reduces 
to  binding  the  solutions  (or  rather  a particular  solution)  to  (lOd), 

In  the  following  section  we  will  look  at  some  ways  of  doing  this. 


21 


C. 


Determination  of 


Our  equation  for 

A.  A 

(H  + U^)  fg  = E 


A / 
B B 


is  of  the  form 


(lla) 


where 


(lib) 


where  we  note  that  fg  is  an  eigenfucntion  of  the  htrmitian  'hamiltonian' 

(H  + U^)  ana  that  H is  the  electronic  hamiltonian  for  isolated  system  B . 
If  we  now  explicitly  assume  that  is  small  compared  to  H (this  of 
course  was  implicit  in  our  derivation  of  U in  the  first  place)  then 
U can  be  regarded  as  being  a small  perturbation  on  H . As  a consequence, 
the  solution,  , to  (11)  which  we  seek  should  resemble  the  electronic 
wavefunction  for  isolated  system  B , fg  , and  this  wavefunction  should 

A 

therefore  provide  us  with  a proper  starting  point  determining  and 

b 


1.  First-Order  Perturbation  Treatment 

Projecting  (lla)  against  fg  and  normalizing  ¥g  to  unity  gives 
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Eg  = <'v'^\H  + ; <'i'g|'t'g>  = 1 , 


(12a) 


whereas  if  S'g  is  a self-consistent  solution  for  isolated  system  B 
its  energy  is 


e;  = = 1 . 


(12b) 


A 

If  we  assume  that  U is  a quite  small  perturbation  to  H then  we  can 
solve  tor  Eg  using  standard  first-order  perturbation  theory  to  obtain 


eJ  = t <.0|uAl,o, 


(13) 


A . 


so  that  our  first-order  perturbative  expression  for  Vg  is  simply 


(14) 


This  of  course  is  equivalent  to  assuming  that  the  wavefunction  for  system 
B remains  virtually  unchanged  when  in  the  presence  of  A . 


In  terms  of  the  one-  and  two-particle  density  matrices,  p°  and  , 
associated  with  Vg  , (14)  becomes 

(15) 

If  Vg  is  a single  Slater  determinant  so  that 
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’’^2’’"i ’’"2^  = l/2[p°(r^  ;rpp°(r2;rp  - p°(ri  •,r2)p°(r2;rp]  (16a) 


P?(r;r')  = 


Z $.  (r)  $,.(r') 
i=l  ^ ^ 


(16b) 


where  {$.}  are  the  n one-electron  spin-orbitals  comprising 

1 ' n 


then  Vg  takes  on  the  form 


’'B.nuc  " '.'2  .2. 


where  the  two-electron  integrals  over  qi2  lire  written  in  the  standard 


<1?||12>  notation. 


If  Hg  cannot  be  written  as  a single  Slater  determinant  but  can 
be  written  as  a linear  combination  of  determinant:  involving  n one-electron 
spatial  orbitals  19°)  then  (12b)  can  always  he  written  as 


Z D^'  <$?|h|i'?>  + z <$?$?|r,„U°$?> 


1 T ' J 


i ,j,k.  I 


^ii 


i kl 

where  are  fixed  coefficients  which  depend  on  the  precise  form 

0 A 

of  Yg  and  on  orbital  overlaps.  Since  U is  of  the  same  form  as  H 

and,  like  H , is  spin  independent  (recall  that  isolated  system  A must 

be  closed-shell)  we  can  immediately  write  our  first-order  approximation 

for  Vg  as 
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B,nuc 


n 

+ L 

1 J 


$?>  + 
J 


i ,J  --k,! 


r«kl  ^0^0 1 A I 

D,j  <»,»j|q,2l 


(19) 


where  all  operators  and  integrals  are  now  explicitly  takent  to  involve 

only  spatial  coordinates.  Therefore,  we  can  determine  <'f'g I U^l 'i'g>  by 

using  the  general  energy  expression  (18)  but  employing  one-  and  two-electron 

integrals  over  the  operators  p^  and  instead  of  the  usual  integrals 

1 

over  h and  r.j2  . 


We  could  continue  this  perturbation  treatment  by  going  on  to  determine 
second-order  and  higher  corrections.  However  we  will  stop  here  and  next 
consider  a self-consistent  approach  instead. 


2.  Self-Consistent  Treatment 

A 

Let  us  now  assume  that  while  U is  a small  perturbation  to  H 
it  is  not  small  enough  to  justify  a simple  first-order  perturbation 

A 

treatment.  That  is,  we  will  now  assume  that  Tg  can  still  be  taken  to 
be  functionally  the  same  as  Yg  but  because  of  the  presence  of  A the 
spatial  orbitals  themselves  distort  away  from  {$9}  to  a meaningful 
extent.  The  problem  is  then  to  determine  these  new  orbitals  . 

Since  Yg  is  an  eigenfunction  of  electronic  motion  satisfying  (11)  this 
can  ba  done  variational ly.  That  is,  we  can  determine  these  new  distorted 

A 

orbitals  by  requirir.g  that  Eg  be  stationary  with  respect  to  changes  in 
these  orbitals.  If  fo*^  simplicity  (but  not  necessity)  we  assume  that 
Yg  can  be  written  as  a hartree-Fock  type  wavefunction  involving  orthonormal 
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spatial  orbitals  (or  orthonortnal  spin-orbitals  for  the  case  of  an 
Unrestricted  Hartree-Fock  wavefunction)  then  (18)  takes  on  the  simple 
Roothaan  form 


Er  = 2 E f.  <«>9|h|$?>  + Z [a..  J^o  0 + b..K  o.o] 

B 1 i'  ' 1 . ij  ij 


where 

J*o.o  = 

1 J 

<<i>9 1 J .0 1 <i>9> 

J 

= «j>9|.)  o| 

J 

<$Vi 

1 J ‘ 

^12  1 J 

K 0 0 = 

«j>9Ik. 

1 ' >{> .'  1 

J 

= 1 
j' 

(j)9>  = 

J 

*0*0  1 
<!{'  .It  . 

1 J' 

^2  J ^ 

(20a) 


(20b) 


(20c) 


and  where  . . ;b . . } are  fixed  coefficients  (for  the  case  of  a multi- 

I I J I J 

configurational  wavefunction  these  coefficients  are  simply  related  to  tne 
variational ly  determined  configurational  coefficients)  which  depend  on 
the  precise  form  of  Y®  . Since  Yg  is  being  taken  to  have  the  same 
form  as  Yg  and  because  of  the  similarities  between  H and  mentioned 
above,  we  can  immediately  write  our  equation  for  Eg  as 


A 

z f .<$•  I h I <!>.>  + 
i=l  ^ ^ ^ 


z 

i ..1 


(21a) 


where  the  operator 


h^  and 


are  given  by 


h^ 


= h + p^ 


(21b) 


Its ii. 
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Therefore,  {'!'.}  and  therefore  Eg  can  be  determined  using  standard 

basis  set  expansion  (LCAO)  SCF  techniques.  The  only  difference  is  that 

instead  of  using  the  usual  one-  and  two-electron  integrals  we  must  use 

A A 

integrals  over  the  operators  h and  g^^  instead.  However,  since 
these  integrals  serve  only  as  imput  this  difference  is  transparent  to 
whatever  available  SCF  procedure  we  employ. 

3.  Configuration  Interaction  Treatment 

For  the  sake  of  completeness  it  should  be  mentioned  that  since 

fl 

Tg  can  be  variationally  determined  there  is  no  need  to  stop  at  the  SCF 
level  and  we  could  solve  for  this  function  as  a Configuration  Interaction 
(Cl)  problem.  The  only  difference  form  a standard  Cl  calculation  is 
that  instead  of  using  the  usual  one-  and  two-electron  integrals  we  must 
use  integrals  over  the  h and  operators  defined  in  (21b).  It 

n 

should  be  pointed  out  however  tha*  any  solution  for  Tg  which  differs 

0 A 

significantly  from  Tg  implies  that  for  that  solution  U can  no  longer 

be  regarded  as  a small  perturbation  to  H and  in  such  a case  the  validity 

n 

of  the  approximations  made  in  our  choice  of  U would  become  subject 
to  question.  This  of  course  also  applies  in  our  SCF  treatment  as  well. 

D.  The  Potentials 

A A 

The  one-  and  two-particle  effective  potentials  E and  W 
appearing  in  our  final  equations  in  section  II  are  hermitian  adiabatic 
approximates  to  the  true  field  theoretic  potentials.  As  we  have  mentioned 


such  potentials  can  be  obtained  in  closed  ab  initio  form  using  many-body 
theory  within  the  framework  of  the  RPA  approximation  and  taking  the 
adiabatic  limits.  However,  while  these  potentials  are  tractable  they  are 
nonetheless  quite  complicated  and  their  use  would  entail  considerable 
computational  effort.  In  view  of  the  perturbative  nature  of  our  theory 
it  is  reasonable  to  expect  that  we  could  use  potentials  having  simpler 
forms.  Such  simpler  potentials  can  be  obtained  by  making  moment  expansions 
of  the  RPA  potentials  and  truncating  these  expansions  in  a physically 
meaningful  manner  [5].  When  this  is  done,  the  resulting  potentials  can 
he  cast  in  forms  which  are  very  similar  to  phenomenologically  derived 
semiempirical  potentials  which  have  been  used  by  other  workers  with 
considerable  success  [6].  Therefore,  it  would  seem  that  the  use  of  com- 
plicated £b  initio  potentials  is  not  warranted  (although  we  do  reserve 
the  option  to  do  so)  and  that  we  can  take  our  potentials  to  have  semi- 
empirical forms  similar  to  those  used  by  Delgarno  and  by  Victor  [6], 
namely 


/dr'  ^ . /dr'  ^(r-.r' ) - W5(k|r| ) 


2|r 


“6  Wg(k|r!)  + (a^  + a^  |r|  + a^kl  )e 


2^  -klri 


(22a) 


and 


(22b) 


Jdr^  dr^  W^(r^,r2;r^,r^)  = - — 2 I *"2 1 ^*^1  ‘^12^ 


|7|J3jr  '13  ■ W4(kkil)W^(k|r2l)P2(cos  Y|2) 
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where  we  have  chosen  our  coordinate  system  to  be  centered  on  A which 
for  simplicity  we  now  take  to  be  an  atom  having  a nuclear  charge  of  z 

I r 


and  where 


and 


Wn(x) 


n2 

P,(x) 

A 

“d 

A 

“q 


(a.) 


= (1-e"^  ) is  a cutoff  function 

= angle  between  vectors  r-j  and  r2 
= legendre  polynomial  of  the  degree 

= dipole  polarizability  of  A 

= approximate  quadrupole  polarizability  of  A ^adjustable) 
= approximately  1/2  r^  where  r^  is  the  'Effective  radius 
of  A (adjustable) 

= adjustable  monopole  parameters 


(r;r') 

HP 


"a  (r')[2-P^^^.]  ^^(r’) 
i=l 


r-r 


/\ 

is  the  static  Hartree-Fock  potential  of  A with  } being  the  n^ 

spatial  Hartree-Fock  orbitals  for  the  electrons  in  A . P„  is  the 

permutation  operator  if  r is  an  electronic  coordinate  whereas  P^  = 0 

I j r 
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if  r is  a nuclear  coordinate.  Note  that  all  quantities  are  now  purely 
A A 

spatial  and  that  T.  and  W are  therefore  explicitly  spin  independent. 

In  (22^)  the  first  term  is  simply  the  potential  due  to  the  nucleus 
of  A and  the  second  term  is  the  static  Hartree-Fock  potential  for  the 

A 

electrons  in  A occupying  the  spatial  orbitals  } . The  next  two 
terms  in  (22a)  are  asymptotically  correct  induced  d pole  and  quadrupole 
polarization  potentials  which  die  off  rapidly  at  short  distances  from  A . 
The  final  term  is  an  induced  monopole  term  which  serves  as  a short  range 
correction  potential.  The  terms  in  (22b)  describe  an  asymptotically 
corre  t dielectric  potential  which  properly  cancels  out  one-particle 
induced  dipole  and  quadrupole  polarizations  of  A due  to  two  particles 
of  the  same  charge  when  these  particles  are  on  opposite  sides  of  A . 
ihat  these  potentials  represent  a significant  simplif'  .ation  ever  the 
ab  initio  potentials  is  clear  in  that  our  two-part i cl e potential  is 
strictly  local  and  the  only  nonlocal  term  in  tne  one-partide  potential 
is  simply  the  usual  riartree-Fock  exchange  potential.  However,  despite 
their  simplicity,  potentials  such  as  these  have  been  used  very  success- 
fully for  a variety  of  problems  in  the  past  and  should  therefore  be  quite 
adequate  for  our  purposes. 

As  it  stands,  (22b)  is  concise  and  to  the  point.  However,  one 
important  property  is  obscured  This  is  that  our  two-particle  potential 
can  be  written  in  terms  of  one-particle  operators  only.  To  see  this 
we  need  only  express  ' ^ angle  y^2 


in  terms  of  the  spherical  angles 
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for  each  vector.  Thus,  by  making  use  of  the  exoansion 

''12>  “ P,(C0S  6,)Pj(cos  ej)  + 2 I fe}f  P”  (MS 

m=l ' 

[cos  m4)^  cos  m<^2  ^in  miti  sin  111412] 

where  ^^(x)  is  an  associated  Legendre  Dolynomial,  we  find 
can  be  rewritten  as 

f A ^ 

Jdr^  dr^  W (r^  ,r2;r]  ,r^)  = - Q^- (r^  )Q^.  (r2) 

where 


Ql(r) 

■W3(klr| )P^ (cos  e) 

Q2(r} 

W3(k| r] ) p| (cos  e)  cos  4’ 

q3(r) 

= 'R  |p|'' 

» 1 

'W3(k|rl ) P^  (cos  e)  sin  4’ 

Q4(r) 

- 

*W4(k!r|)  P2(cor  e) 

Q5(r) 

= ’/3^  1 

rl  ^W,(klrl)  p3,(cos  e)  cos  41 

• u 

Q6(r’) 

' '/3^  1 

r!'^W^(klrl)  P2(cos  0)  sin  4> 

Q7(r) 

= ^nz^R 

Ir]  ^W^(klrl)  P2(cos  0)  cos  241 

%{r) 

- 1/12^ 

lr|  \,(k|r!)  P2(cos  e)  sin  24> 

^)p'J(cos  62) 
(24) 

that  (22b) 

(25a) 


(25b) 


31 


1. 


4. 


5. 

6. 


1).  Schneider,  H.  S.  Taylor  and  R. 

G.  Csanak  and  H.  S.  Taylor,  Phys.  Rev.  A6,  1843  (197  ). 


B.  S.  Yarlagadda,  G.  Csanak,  H.  S.  TAylor,  R.  Schneider  ^ Yaris, 
Phys.  Rev.  A7,  146  (1973);  L D.  Thomas,  G.  Csanak,  H.  S.  Taylor  an 
B.  S.  Yarlagadda,  J.  Phys.  B7,  1719  (1974). 


B.  Schneider,  H.  S.  Taylor,  R.  Yaris,  and  B.  S.  Yarlagadda,  Chem. 
Phys.  Letters  381  (1973). 


G.  Csanak  and  H.  S.  Taylor,  Phys. 

G.  Csanak  and  H.  S.  Taylor,  J.  Phys.  B6,  2055  (1973). 


S.  W.  Wang,  H.  S.  Taylor  and  R.  Yaris,  Chem.  Phys.  1_4,  53  (1976). 


A.  Dalqarno,  C.  Bottcher  and  G.  A.  Victor,  Chem.  Phys.  Lett.  , 

265  (1970),  C.  Bottcher  and  A.  Dalgarno,  Proc.  R. 

187  (1974);  D.  K.  Watson,  C.  J.  Cerjan,  S.  Gubermen  and  A.  Dalgarno, 

submitted  to  Chem.  Phys.  Lett.  1977.  nq7Bl 

C.  Laughlin  and  G.  A.  Victor,  Atomic  Physics  3,  247 


